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Inspired by the recent theoretical suggestion that the random-bond S' = 1/2 antiferromagnetic 
Heisenberg model on the triangular and the kagome lattices might exhibit a randomness-induced 
quantum spin liquid (QSL) behavior when the strength of the randomness exceeds a critical value, 
and that this “random-singlet state” might be relevant to the QSL behaviors experimentally ob¬ 
served in triangular organic salts k—(ET) 2 Cu 2 (CN) 3 and EtMe 3 Sb[Pd(dmit) 2]2 and in kagome her- 
bertsmithite ZnCu 3 (OH) 6 Cl 2 , we further investigate the nature of the static and the dynamical spin 
correlations of these models. We compute the static and the dynamical spin structure factors, S(q) 
and 5(q,aj), by means of an exact diagonalization method. In both triangular and kagome mod¬ 
els, the computed ^(q, cu) in the random-singlet state depends on the wavevector q only weakly, 
robustly exhibiting gapless behaviors accompnied by the broad distribution extending to higher 
energy uj. Especially in the strongly random kagome model, S(q,a;) hardly depends on q, and 
exhibits an almost flat distribution for a wide range of uj, together with a, uj — 0 peak. These 
features agree semi-quantitatively with the recent neutron-scattering data on a single-crystal her- 
bertsmithite. Furthermore, the computed magnetization curve agrees almost quantitatively with 
the experimental one recently measured on a single-crystal herbertsmithite. These results suggest 
that the QSL state observed in herbersmithite might indeed be the randomness-induced QSL state, 
i.e., the random-singlet state. 


I. INTRODUCTION 

Geometrically frustrated magnets have attracted a 
long-standing and ongoing attention in the field of con¬ 
densed matter physics because they often give rise to 
a variety of nontrivial thermodynamic states. In par¬ 
ticular, the quantum spin liquid (QSL) state having no 
magnetic long-range order (LRO) has been extensively 
investigated both theoretically and experimentally. An¬ 
derson proposed the resonating valence bond (RVB) state 
as a possbile ground state of the S' = 1/2 Heisenberg an- 
tiferromagnet on the triangular lattice^. According to 
many subsequent theoretical studies, however, it is now 
widely believed that the ground state of a simple S = 1/2 
antiferromagnetic (AF) Heisenberg model on the trian¬ 
gular lattice with the nearest-neighbor bilinear interac¬ 
tion exhibits a Neel LRO with the three-sublattice 120° 
structure^Ti^. 

By contrast, the ground state of the S = 1/2 Heisen¬ 
berg model on the kagome lattice, which has stronger 
frustration than the triangular lattice, is quite likely to 
be non-magnetic. Namely, the absence of the magnetic 
LRO in the S = 1/2 kagome AF Heisenberg model has 
now been established from various numerical studies^^— , 
although the precise nature of its ground state still re¬ 
mains controversial. Various candidates, including the 
gapped Z 2 spin-liquidiid^i^i^^, the gapless U(l) spin 
liquidi^i^^— , the chiral spin liquid^^ and the valence 
bond crystaliii^Si^ etc, have been proposed. 

Along with these theoretical studies, experimen¬ 
tal quest for the QSL has also been persued. As 
a result, several promising candidate materials were 
recently reported. In the triangular system, cer¬ 


tain organic salts such as k— (ET) 2 Cu 2 (CN) 3 ^“— and 
EtMe 3 Sb[Pd(dmit) 2 ] 2 ^— were reported to exhibit a 
QSL-like behavior, i.e., exhibit no magnetic LRO down 
to a very low temperature. These organic materials 
show gapless (or nearly gapless) behaviors at low tem¬ 
peratures. For example, the specific hea t^^d*^ and the 
thermal conductivit} ^^^^ exhibit the behavior linear in 
the absolute temperature T. Obviously, in understanding 
the QSL-like properties of these triangular organic salts, 
some extension beyond the simplest nearest-neighbor 
Heisenberg model is needed^^— . 

In the kagome system, herbertsmithite 
ZnCu 3 (OH) 6 Cl 2 was reported to be a candidate of 
QSL, again showing gapless behaviors in various phys¬ 
ical quantities'^— . The true nature of these gapless 
QSL candidates experimentally observed both in the 
triangular and the kagome lattice AFs has still remained 
obscure in spite of much theoretical and experimental 
efforts devoted to the issue. 

For years, major part of theories have tried to elucidate 
the experimentally observed QSL behaviors as properties 
of a clean and regular system. By contrast, it was pro¬ 
posed recently in Refs— and^ that the quenched bond¬ 
randomness, together with the geometrical frustration 
effect, might be essential in stabilizaing the experimen¬ 
tally observed gapless QSL state both in the triangular 
organic salts and in the kagome herbersmithite. Such 
a randomness-induced QSL state is called the “random- 
singlet state”— or the “valence bond glass state”— 
where the spin singlet is formed in a spatially ran¬ 
dom manner. Possible importance of the quenched site- 
randomness in kagome herbertsmithite was also pointed 
out in Ref—. 
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The origin, or even the existence of the quenched ran¬ 
domness in triangualr organic salts and kagome herber- 
smithite is a non-trivial matter, especially in view of the 
fact that the theory requires a considerable amount of 
randomness, not an infinitesimal one, to realize the QSL- 
like behavio r^^'^’^ . 

In the case of the triangular organic salts, it was sug¬ 
gested that the randomness for the spin degrees of free¬ 
dom was self-generated at low temperatures via the ran¬ 
dom freezing of the electric-polarization degrees of free¬ 
dom inherent to these organic salts consisting of molecu¬ 
lar dimers^. In fact, the measured ac dielectric constant 
of these organic salts exhibited a glassy response even 
on a macroscopic time scale of kHz in the temperature 
region where the QSL behavior was observed in the spin 
degree of freedom^^. 

In the case of kagome herbertsmithite, the quenched 
randomness comes from the random substitution of non¬ 
magnetic Zn^+ by magnetic Cu^+, located on the tri¬ 
angular layer adjacent to the kagome laye r^^i^^ . Note 
that the kagome layers in herbertsmithite are separated 
by [Zn(OH)6]"‘“ octahedral units whose Zn^+ constitutes 
the triangular layer. It was reported that about 15% of 
Zn^+ on the triangular layer was randomly substituted by 
Cu^+, keeping the kagome layer intaclj^. Since Cu^'*' is 
a Jahn-Teller ion, such a random substituion would lead 
to a random Jahn-Teller distortion of the [Cu(OH)6]^“ 
octahedra, leading to the random modification of the ex¬ 
change path, and subsequently the exchange strength, 
connecting the Cu^’*' on the kagome layer—. 

Indeed, the numerical results on a simplihed random 
S' = 1/2 AF Heisenberg model on the triangular and 
the kagome lattices appear to reproduce many of ex¬ 
perimentally observed feature o^^’^^ of various thermody¬ 
namic quantities, including the T-lienar low-temperature 
specific hea t^^d^’^^ , the gapless magnetic susceptibility 
occasionally accompanying a Curie-like tail^, and the 
gapless temperature dependence of the NMR relaxation 
rate 

In view of this apparent success of the random model 
in reproducing the experimentally observed QSL-like be¬ 
haviors, it would be desirable to further investigate the 
nature of the static and the dynamical spin correlations of 
the randomness-induced QSL state, the random singlet 
state. For this purpose, in the present paper we com¬ 
pute by means of an exact diagonalization technique the 
static and the dynamical spin structure factors which are 
experimentally accessible via, e.^., the elastic and the in¬ 
elastic neutron scattering measurements. Such a compar¬ 
ison between theory and experiment might give further 
information in examining the validity of the randomness- 
induced QSL picture of the experimentally observed QSL 
states. In order to clarify the effect of frustration, we per¬ 
form comparative calculations also on the random-bond 
S' = 1/2 AF Heisenberg model on the square lattice. 

Our numerical results corroborates the previous obser¬ 
vation that both the triangular and the kagome models 
exhibit the randomness-induced QSL-like behavior when 


the randomness exceeds a critical value^i^. Meanwhile, 
the unfrustrated square model persistently exhibits the 
AF LRO up to the maximal randomness without showing 
the QSL-like behavior. The result hilights an important 
role of frustration, along with the randomness and the 
quantum fluctuation, in stabilizing the random-singlet 
state. 

The random-singlet states in the triangular and in the 
kagome models have some mutual similarities, but also 
some differences. In the triangular case, the random- 
singlet state keeps a certain amount of AF short-range 
order even at the maximal randomness. While the dy¬ 
namical structure factor S{q,uj) of the triangular model 
exhibits a signature of the AF LRO and the magnon 
excitation in the regular and weakly random cases, it 
exhibits, in the strongly random case corresponding to 
the random-singlet state, a gapless behavir accompa¬ 
nied by a broad w-distribution extending to higher en¬ 
ergy, which is dependen on the wavevector q only weakly. 
In the kagome case, by contrast, a signature of the AF 
LRO or the magnon excitation is hardly discernible ei¬ 
ther in the regular or in the random case. Peaky fea¬ 
tures of S'(q, uj) still retained in the regular case give way 
to gapless behaviors in the strongly random case cor¬ 
responding to the random singlet state. Such features 
of 5(q, w) are accompanied by an almost flat distribu¬ 
tion in a wide range of lo and by an a; = 0 peak, which 
hardly dependes on the wavevector q. Indeed, these fea¬ 
tures of the computed S'(q, w) of the strongly random 
kagome model are compared quite faborably with the re¬ 
cent inelastic-neutron scattering data on a single-crystal 
kagome herbersmithite^. 

This paper is organized as follows. In Sec. H, we 
present our model and the details of our numerical cal¬ 
culation. In Sec. HI and IV, we show the results of our 
numerical calculations in the case of the triangular and 
the kagome models, respectively. Section V is devoted 
to summary and discussion. For comparison, we also 
present the corresponding results for the unfrustrated 
square-lattice model in Appendix A. The detailed infor¬ 
mation about the shape of finite-size lattices employed 
in our exact-diagonalization calculation is given in Ap¬ 
pendix B. 

II. THE MODEL AND THE METHOD 

Our model is the AF random-bond S' = 1/2 quantum 
Heisenberg model on the triangular and the kagome lat¬ 
tices, whose Hamiltonian is given by 

n = J2 ■ s. (1) 

id 

where = (Sf, Sf, Sf) is a spin-1/2 operator at the *-th 
site on the lattice, while Jij > 0 is the random nearest- 
neighbor AF coupling obeying the bond-independent uni¬ 
form distribution between [(1 — A) J,(l -|- A) J] with the 
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Figure 1. (Color online) The temperature dependence of 
the rescaled sublattice magnetization per spin fhs of the 
triangular-lattice Heisenberg antiferromagnet for the random¬ 
ness A = 0,0.7 and 1.0 and for the sizes = 12 and 18. The 
definition of ms is given in the text. 


mean J. The parameter A represents the extent of the 
randomness: A = 0 and A = 1 cases correspond to the 
regular and maximally random cases, respectively. 

According to Refs.— and^, the random-singlet ground 
state is realized when the randomness is stronger than a 
critical value Ac. In the triangular model, Ac is esti¬ 
mated to be ~ 0.6 where Ac separates the AF phase 
and the random-singlet phased, while in the kagome 
model it is estimated to be Ac ~ 0.4 where Ac separates 
the the randomness-irrelevant QSL phase {e.g., the Z 2 
spin-liquid phase) and the randomness-relevant random- 
singlet phased. 

In the present paper, we employ the exact diagonaliza- 
tion (ED) method in computing various physical quanti¬ 
ties. The ED method is precise and is applicable even to 
systems with frustration, while it has a disadvantage of 
being limited to very small system sizes. In our compu¬ 
tation, the total number of spins N is N = 12,18, 24,30 
for T = 0 and N = 12,18 for T > 0, periodic boundary 
conditions being employed. Sample average is taken over 
100 {N = 12,18, 24) and 50 {N = 30) independent bond 
realization in the T = 0 calculation, while 100 {N = 12) 
and 30 {N = 18) in the T > 0 calculation in both cases of 
the triangular and the kagome models. The shape of the 
lattice is illustrated in Fig. 13 and Fig. 14 in Appendix B 
for the triangular and the kagome lattices, respectively. 
In what follows, the energy and temperature are normal¬ 
ized in units of J. 


III. RESULTS I: THE TRIANGULAR LATTICE 

In this section, we study the ground-state and the 
finite-temperature properties of the random-bond S = 
1/2 AF Heisenberg model on the triangular lattice. We 
begin with the sublattice magnetization associated with 
the three-sublattice 120° structure. The sublattice mag¬ 


netization per spin nis may be defined by 





1 

(A/6)(A/6 + l) 



12 ^ 

“ N{N + 6) ^ 


i,j^a 


( 2 ) 


where a = 1,2,3 denotes three triangular sublattices 
of the original triangular lattice, and the sum i £ a 
{i,j € a) is taken over all site i (i and j) belonging to 
the a-th sublattice. The symbol (• • •) denotes the ther¬ 
mal average, and [• • • ] the average over the bond disorder. 
Note that m-s is normalized to give unity for the classical, 
perfectly ordered 120° structure. 

We also define the rescaled sublattice magnetization 
per spin mg by subtracting the auto-correlated part as. 




12N 3 

N(N + 6) ^ 4’ 


( 3 ) 


The subtraction of the auto-correlated part is made so 
that the sublattice magnetization tends to zero at T —^ 00 
even for finite N. In the thermodynamic limit N —>■ 00 , 
fhs approaches rus at any temperature. 

In Fig. 1, we show the temperature dependence of 
the rescaled sublattice magnetization per spin, mg, for 
N = 12 and 18 for several values of the randomness A. In 
the regular and weakly random cases, fhs increases mono- 
tonically with decreasing the temperature, and eventu¬ 
ally saturates. This increase of rUg at finite temperature 
is associated with the growth of the AF short-range or¬ 
der. When the randomness gets stronger, fhs is supressed 
as a result of the supression of the AF short-range order 
due to the randomness. 

An interesting observation is that, for the strongest 
randomness A = 1, fhs tends to decrease weakly with de¬ 
creasing the temperature in the lower temperature range 
of T < 0.2, exhibiting some sort of crossover. This su¬ 
pression is weak, but tends to be more eminent for larger 
sizes. Such supression of the AF short-range order is com¬ 
patible with the formation of the random-singlet-type 
state at T < 0.2. Similar supression was also reported 
in Refi^ in the temperature dependence of l/Ti in the 
similar temperature and size range. There, Z 2 vortesi^^ 
was invoked to be a possible candidate of this crossover¬ 
like anomaly, since the supression was observed only for 
lattices larger than N = 18. Further study is required to 
clarify the nature of this crossover phenomenon, which 
is observed only for modestly large systems with strong 
randomness. 

In order to get further information about the static 
spin correlations of the model, we investigate the static 
spin structure factor S'(q) defined by 


^(q) 


1 

A 




( 4 ) 
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Figure 2. (Color online) Intensity plots of the static spin 
structure factors ^(q) of the triangular-lattice Heisenberg an- 
tiferromagnet in the wavevector {q^, qy) plane for several val¬ 
ues of the randomness A = 0 (a), 0.3 (b), 0.6 (c), and 1.0 (d). 
The lattice constant a = 1 is the nearest-neighbor distance of 
the triangular lattice. The system size is N — 30. The solid 
black line depicts the first Brillouin zone of triangular lattice. 
Each small hexagon corresponds to the resolution unit, and 
the center of each distorted hexagon is the wavevector point 
we can treat in A = 30 system. The black point in Fig. (a) 
represents the K point, while the orange point the M point: 
see the text for details. 


where (• • •) means the ground-state expectation, q is the 
wavevector, and Rj is the position vector at the site j. 
The computed static spin structure factor at T = 0 is 
shown in Fig. 2 as an intensity plot in the wavevector 
{Qx, Qy) plane for several values of the randomness, A = 0 
(a), 0.3 (b), 0.6 (c), and 1.0 (d). The system size is 
A = 30. 

Although the resolution is rather low due to the small 
system size, several characteristic features are clearly dis¬ 
cernible from the figure. In the regular case A = 0, emi¬ 
nent peaks corresponding to the AF LRO are observed at 
the so-called K points, q = {qx,Qy) = (±27r/3, ±27r/-\/3) 
and (±47r/3, 0), where the length unit is taken as the 
nearest-neighbor distnace of the triangular lattice. As 
the randomness gets stronger, the peak at the K points 
is gradually suppressed. The AF LRO is expected to 
vanish for A > Ac ~ 0.6^. Such a phase-transition-like 
sharp change of behavior around Ac ^ 0.6, however, is 
not necessarily clear here, presumably due to the small 
system size of A = 30. In fact, a rounded peak per¬ 
sists at the K point even at the maximal randomness of 
A = 1, suggesting the persistance of the AF short-range 
order in the random-singlet state. Thus, in the triangu¬ 
lar case, the random-singlet state coexists with the AF 
short-range order. 


In order to get further information about the dynam¬ 
ical spin correlations of the model, we next investigate 
the dynamical spin structure factor defined by 

S{<i,w) = - (£;„ - Eo))] , 

n 

(5) 


where is the z-component of the Fourier transform of 
the spin operator, ipn is the eigenfunction of the Hamil¬ 
tonian (I) whose eigenvalue is A„, and V'o is the ground- 
state eigenfunction with the eigenvalue Eq. By using the 
continued fraction method^^, it may be rewritten as 


S'(q, w) = — lim —Im 

TT 


- lim — Im 

ly-s-O TT 




I 


uj Eq — V.irj 




uj -\- Eq-\-ir] — aq— 




where ai and are the diagonal and sub-diagonal el¬ 
ements of the tri-diagonal Hamiltonian obtained by the 
Lanczos method. In implementing the continued-fraction 
expansion, we performed at least 1000 iterations. A small 
but finite r], rj = 0.02, is used. 

The oj-dependence of 5'(q, w) computed at the K point 
is shown in Fig.3, while the one at the M point, which is 
located at the midpoints of the 1st Brillouin zone (BZ) 
edges, i.e., q = (0, ±27r/-\/3) and (±7r, ±7r/-\/3), is in 
Fig.4. The randomness A is taken to be A = 0 (a), 0.3 
(b), 0.6 (c), and I.O (d). Note the difference in the ordi¬ 
nate scale between Figs. 3 and 4, and between Figs, (a) 
and (b)-(d). 

Reflecting the AF short-range order of the triangular 
model, the S{q,u!) intensity tends to be larger at the K 
point (Fig.3) than at the M point (Fig.4) irrespective of 
the A-value. The contrast between the K and the M 
points, however, tends to be milder for larger A. 

In the regular case A = 0, the dominant peak observed 
in 5'(q, w) at the K point (Fig.3(a)) is a single magnon 
excitation. Indeed, the peak location in uj tends to zero 
when N is increased toward the thermodynamic limit. 
Similar behavior is observed also in the weakly random 
case of A = 0.3: see the insets of Figs.3(a) and (b). As 
shown in Fig.4(a), several peaks are observed at the M 
point. According to the recent spin-wave analysis for the 
regular model, S'(q, w) at the M point exhibits a leading 
peak at w ~ 0.8^1. Indeed, if one looks at our data of 
A = 0 for A = 24 and 30, the dominant peak appears at 
a nearby position w ~ I. 

As the randomness becomes stronger getting into the 
random-singlet phase, 5'(q, w) exhibit less peaky behav¬ 
ior both at the K and the M points. As can be seen from 
Figs.4(c) and (d), it exhibits a broad distribution extend¬ 
ing to higher w, with a finite intensity gradually growing 
at w = 0. This demonstrates that the random-signlet 
state of the triangular model is indeed magnetically gap¬ 
less both at the K and the M points. The observed gap¬ 
less feature of ^(q, ui) is fully consistent with the gapless 
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(C) K(A=0.6) (d) K(A=1.0) 



CO CO 


(a) M(A=0.0) (b) M(A=0.3) 




Figure 3. (Color online) The cj-dependence of the dynamical 
structure factors of the triangular-lattice Heisenberg antifer- 
romagnet taken at the K point, q = (±27r/3, ±27r/\/3) and 
(±47r/3,0), for several values of the randomness A = 0 (a), 
0.3 (b), 0.6 (c), and 1.0 (d). The critical randomness sepa¬ 
rating the AF state and the random-singlet state is Ac — 0.6. 
Note the difference in the ordinate scale between Figs, (a) and 
(b)-(d). In the inset of Figs.(a) and (b), the cd value of the 
dominant peak is plotted versus the Inverse system size 1/A. 


Figure 4. (Color online) The cj-dependence of the dynam¬ 
ical structure factors of the triangular-lattice Heisenberg an- 
tiferromagnet taken at the M point, q = (0, ±27r/%/3) and 
(±7r, ±7r/%/3), for several values of the randomness A = 0 
(a), 0.3 (b), 0.6 (c), and 1.0 (d). The critical randomness 
separating the AF state and the random-singlet state is es¬ 
timated to be Ac — 0.6. Note the difference in the ordinate 
scale between Figs.(a) and (b)-(d), and between this figure 
and Fig.3. 


behavior observed in several other observables in Refi^. 
While the intensity tends to be larger at the K point 
than at the M point, the difference tends to be smaller 
for larger A as in the case of the static spin structure 
factor. 

In the random-singlet state, ^(q, w) exhibits a tail in 
Lj in the higher-w range. The asymptotic w-dependence 
of this tail is found to be exponential ~ exp[— w/wq] with 
a characteristic energy scale ujq . Estimates of ujq yields a 
value around 2~2.5. 


IV. RESULTS II: THE KAGOME LATTICE 

Next, we deal with the random-bond S' = 1/2 AF 
Heisenberg model on the kagome lattice whose Hamil¬ 
tonian is given by eq. (1). We first examine the temper¬ 
ature dependence of the two representative types of the 
rescaled sublattice magnetization per spin, mg, each as¬ 
sociated with the q = 0 and the x structures. As 
in the triangular case, it is normalized to give unity for 
the classical, perfectly ordered g = 0 or -x/S x y/S struc¬ 
ture, while the auto-correlated part is subtracted. The 
temperature dependence of the computed thg is shown in 
Fig.5 for N = 18 for several values of the randomness A. 

Numerically, it has been established that the regular 
model exhibits neither the q = 0 nor the X LRO 
even at T = 0^. These AF orders are not realized in the 
random model, either—. Yet, the temperature and the 
size dependence of mg is expected to provide us useful 


information about the associated AF short-range order. 

For the regular model, while the spin-wave-type 1/S 
expansion suggested the possible dominance of the \/3 x 
•\/3 orderi22i2S, the recent numerical results from the 
exact-diagonalizationH and the DMRG— calculations 
suggested the dominance of the q = 0 state. Though our 
present data for N = 30 apparently suggest the domi¬ 
nance of the X short-range order, our maximum 
size A = 30 is smaller than the ones of Ref— andi^, i.e., 
N = 36 and N = 108, repsectively, and might be subject 
to stronger finite-site effect. In any case, as the random¬ 
ness gets stronger, the difference between the q = 0 order 
and the x '/S order tends to be negligible. 

To get further information about the spin correlations 
of the model, we compute the T = 0 static spin structure 
factor, and the results for N = 30 are shown in Fig.6 as 
an intensity plot in the (qx, qy)-plane for the randomness 
of A = 0 (a), 0.4 (b), 0.7 (c), and 1.0 (d). The length 
unit here is taken to be the nearest-neighbor distance of 
the original kagome lattice. 

In the regular case shown in Fig.6(a), the intensity 
appears primarily along the zone boundary of the ex¬ 
tended BZ, consistently with the behavior reported in 
Refs— We note that the recent numerical studies for 
larger system of A = 36 reveal additional small peaks 
at the wavevectors corresponding to the q = 0 state, 
q = (0, ±27r/-\/3) and (±7r, ±7r/-\/3), which are located 
along the zone boundary of the extended al¬ 

though this peak is not clearly discernible in our present 
data of A = 30. 

As can be seen from Figs.6(a)-(d), qualitative features 
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(a) A=0.0 (b) A=0.3 




Figure 5. (Color online) The temperature dependence of the 
rescaled sublattice magnetization per spin, rhs, associated 
with the q = 0 order and the %/3 x %/3 order of the kagome- 
lattice Heisenberg antiferromagnet for several values of the 
randomness A = 0 (a), 0.4 (b), 0.7 (c), and 1.0 (d). 


of S'(q) do not change much even when the randomness 
is introduced up to A = 1, except that the ridge-like 
intensity along the extended BZ boundary is somewhat 
broadend. 

We also compute the dynamical spin structure factor 
S'(q, w) at the two representative wavevectors, i.e., the 
so-called F point, q = (0, ±27r/-\/3) and (±7r, ±7r/-\/3), 
and the M point, q = (0,±7r/-\/k) and (±7r/2, ±7r/2-\/3). 
The computed w-dependence of ^(q, w) is shown for sev¬ 
eral values of the randomness A in Fig.7 for the F point, 
and in Fig.8 for the M point. The w-dependence of 
the dynamical spin structure factor shows some differ¬ 
ences between in the randomness-irrelevant QSL state at 
A < Ac and in the randomness-relevant QSL state at 
A > Ac. In contrast to the case of the triangular model, 
S'(q, w) exhibits a rather broad distribution even in the 
regular case, yet with some pronounced peaks remaining 
both at the F- and the M points as shown in Figs.7(a) 
and 8(a). Although our data of A < 30 are still subject 
to considerable finite-size effects, Lauchli et al. reported 
that S'(q, w) for iV = 24 and for larger A = 36 came 
close, suggesting that 5'(q, w) for A = 30 were not far 
from that of the bulk system. The issue of a small non¬ 
zero gap exists or not in the regular systems is beyond 
the capability of our present calculation. 

For larger A > Ac ~ 0.4 corresponding to the random- 
singlet state, by contrast, ^(q, w) exhibits both at the F 
and M points a nealy flat distribution in the wide w-range 
of, say, w ;< 1.5 and extends to higher w, with a nonzero 
intensity growing at w = 0. This domonstrates the gap¬ 
less nature of excitations in the random-singlet state both 
for the F- and M points. In the case of the strongest ran¬ 
domness A = 1.0, even an w = 0 peak appears at both 
F and M points. As can be seen from Fig.7(d) and 8(d), 


(a)A=0.0 


(b)A=0.4 




Figure 6. (Color online) Intensity plots of the static spin 
structure factors S'(q) of the kagome-lattice Heisenberg anti¬ 
ferromagnet in the wavevector (qx , qy ) plane for several values 
of the randomness A = 0 (a), 0.4 (b), 0.7 (c), and 1.0 (d). 
The lattice constant a = 1 is the nearest-neighbor distance 
of the kagome lattice. The system size is A = 30. The solid 
black line depicts the zone boundary of the extended BZ. The 
black square in Fig.(a) represents the F point, while the blue 
triangle represents the M point. 


the overall behavior of iS'(q, w) are very much similar at 
the F and the M points in the random-singlet state. 

In the random-singlet state, ^(q, w) exhibits a tail in 
Lo in the higher-w range both at the F and the M points. 
The asymptotic w-dependence of this tail is found to be 
exponential, ~ exp[—w/wo], with a characteristic energy 
scale wq. In the kagome model, such an exponential tail 
is also realized even in the regular system. Estimates of 
Wo yields a value around 1.5^2, slightly increasing with 
increasing the randomess A. 

The behavior of 5'(q, w) in the random-singlet state 
of the kagome model is rather similar to the one of the 
triangular model shown in Figs. 3 and 4. One difference 
might be that 5'(q, w) in the low-energy region is even 
more flatter in the kagome model, and the w = 0 peak is 
more eminent. 

Anyway, the gapless behavior with an almost flat dis¬ 
tribution extending to higher w irrespective of the q-value 
is a common feature of the observed dynamical structure 
factor of both the triangular and the kagome models for 
larger A, and might be regarded as a characteristic of the 
random-singlet state. 

In order to make comparison with the recent in¬ 
elastic neutron-scattering data on a single-crystal of 
herbertsmithite^, we provide in the insets of Figs. 7 and 
8 their low-energy part of the experimental data, which 
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Figure 7. (Color online) The tj-dependence of the dynamical 
spin structure factor ^(q, oj) of the kagome-lattice Heisenberg 
antiferromagnet taken at the F point, q = (0, ±27r/'\/3) and 
(±7r, Ttt/v^), for several values of the randomness A = 0 (a), 
0.3 (b), 0.6 (c), and 1.0 (d). The critical randomness sepa¬ 
rating the randomness-irrelevant QSL stae and the random- 
singlet state is estimated to be Ac — 0.4. The insets are 
magnified views of the low-m region in units of meV, where 
J = 17 meV is assumed with herbertsmithite in mind. Note 
the difference in the ordinate scale between Figs.(a) and (b)- 
(d). 


corresponds to the range indicated by the arrow in the 
main panel. Since the experimental data now available 
are limited to the low-energy region of w < lOmeV, we 
show in the inset exactly this range, z.e., u! < 10 meV, 
by using an experimental estimate of J ~ 17meV. This 
estimate of the exchange coupling J of hermertsmithite 
is made via neutron-scattering measurements on powder 
samples^. On comparison with the corresponding ex¬ 
perimental data (Fig.2 of Refi^), one finds a nice agree¬ 
ment especially for the strong randomness A = 1 shown 
in Fig.(d), including the features of (1) a plateau-like be¬ 
havior observed for 2 < w < 10 meV, (2) an a; = 0 
peak observed at w < 1.5 meV without any gap, and (3) 
the insensitivity of the overall behavior of S'(q, w) on the 
wavevector q. We emphasize that the data of the random 
model shown in Fig.(d) appears to resemble the experi¬ 
mental one much more than that of the non-random sys¬ 
tem of Fig.(a). This observation certainly lends support 
to the view that the randomness is playing a signihcant 
role in the QSL state of herbertsmithite. 

Very good agreement with experiment is also found for 
the magnetization curve. In Fig. 9, we show the ground- 
state magnetization curve for the randomness A = 0.7 
for TV = 18 and 24, where the magnetization M is nor¬ 
malized per formula unit of herbertsmithite, i.e., the sat¬ 
uration value taken to be 3, while an applied field is given 
in units of tesla with assuming the experimental g-factor 
{g ~ 2.2) and J-value (J ~ 200K) of herbertsmithite^. 
Then, the computed magnetization curve of Fig.9 indeed 




Figure 8. (Color online) The tj-dependence of the dynamical 
spin structure factor 5'(q, ca) of the kagome-lattice Heisen¬ 
berg antiferromagnet taken at the M point, q = (0,±7r/%/3) 
and (±7r/2, ±7r/2%/3), for several values of the randomness 
A = 0 (a), 0.3 (b), 0.6 (c), and 1.0 (d). The critical random¬ 
ness separating the randomness-irrelevant QSL stae and the 
random-singlet state is estimated to be Ac — 0.4. The in¬ 
sets are magnified views of the low-ca region in units of meV, 
where J = 17 meV is assumed with herbertsmithite in mind. 
Note the difference in the ordinate scale between Figs, (a) and 
(b)-(d). 



H[T] 

Figure 9. (Color online) The T = 0 magnetization curve 
for the randomness of A = 0.7 for A = 18 and 24, in the 
wider field range (inset) and in the lower field range (main 
panel) corresponding to the region indicated by the arrow in 
the inset. With herbertsmithite in mind, the magnetization 
M is normalized per formula unit of herbertsmithite, i.e., 
the saturation value taken to be 3, while an applied held is 
given in units of tesla with assuming the experimental y-factor 
{g ~ 2.2) and J-value (J ~ 200K) of herbertsmithite^. 


exhibits an almost quantitative agreement with the re¬ 
cent experimental data given in Fig. 2 of Refi^ in the 
same units. 

Two features might be noticed. One is the absence of 
any plateau-like anomaly in the magnetization curve. As 
shown in the inset, it exhibits a near linear behavior in 
































































an entire region of H up to the near saturation, except 
for some ‘wavy’ behavior occurring at intermediate fields 
(which becomes less visible for the stronger randomness 
of A = 1—). As reported Refi^, we have found that the 
plateau-like anomaly tends to go away for A > Ac, yield¬ 
ing a near linear behavior. The other notable feature of 
the magnetization curve might be an upper-convex gap¬ 
less behavior observed at weaker fields near H = 0. This 
enhanced behavior of the low-field magnetization is con¬ 
sistent with the Curie-like behavior of the susceptibility 
observed for stronger randomness^, and is likely to be 
borne by the ‘free’ or ‘almost free’ spins inevitably gen¬ 
erated in the random-single state. 


V. SUMMARY AND DISCUSSION 

We studied by means of an ED method the nature 
of spin correlations of the random-bond S' = 1/2 AF 
Heisenberg models on the triangular and the kagome lat¬ 
tices via the static and the dynamical spin structure fac¬ 
tors. To hilight the possible importance of frustration, 
we also made a comparative calculation for the unfrus¬ 
trated random-bond S = 1/2 AF Heisenberg model on 
the square lattice. 

Both the triangular and the kagome models exhibit the 
randomness-induced QSL behavior when the randomness 
exceeds a critical value as observed in previous studies, 
while the unfrustrated square model persistently exhibits 
the AF LRO up to the maximal randomness without 
showing the QSL behavior. This demonstrates that the 
frustration is certainly playing a role in stabilizing the 
random-singlet state. 

The random-singlet states in the triangular and the 
kagome models have some similarities, but also some dif¬ 
ferences. The random singlet state of the both models 
exhibit gapless behaviors, dependent on the wavevector 
q only weakly, while the dynamical spin structure fac¬ 
tor S'(q, w) exhibit a broad distribution in uj extending 
to higher lj with an exponentail tail. Especially in the 
strongly random kagome model, /S'(q, w) hardly depends 
on q and exhibits an almost flat distribution for a wide 
range of lo with a w = 0 peak. 

As discussed in §IV, our results for the dynamical 
spin structure factor for the strongly random kagome 
model compares quite favorably with the recent inelastic 
neutron-scattering data on the kagome herbersmithite, 
including (1) a plateau-like behavior observed for 2 < 
O' < 10 meV, (2) an a; = 0 peak observed w < 1.5 
meV without any gap, and (3) the insensitivity on the 
wavevector q. Since the present experimental data are 
limited to the low energy range w ;< 10 meV, it would be 
interesting to perform further experiments in the higher 
energy range of w > 10 meV to make further compari¬ 
son. In addition, the computed magnetization curve is 
found to exhibit a good, almost quantitative agreement 
with the recent experimental data on herbertsmithite. 

In making a truly quantitative comparison with the 


experimental data on herbertsmithite, ones needs to ex¬ 
amine several effects not considered in the present model, 
i.e., the effect of the Dzyaloshinskii-Moriya interaction 
and the triangular layer between the kagome layers, etc. 
How the results depend or do not depend on the partic¬ 
ular form of the randomness needs further clarification. 
In particular, the possible effect of the dilution-type ran¬ 
domness pointed out experimentally in RefSi ^°i^^i^^ and 
studied theoretically in Ref<2^ might be incorporated. 

An important open theoretical question in the kagome 
model might be the distinction and the true relation be¬ 
tween the observed two types of QSL-like states, i.e., 
the randomness-irrelevant QSL state realized at A < Ac 
and the randomness-relevant QSL state (random-singlet 
state) realized at A > Ac. One intuitive view of this 
transition (or crossover) occurring at A ~ Ac might be 
that singlet states tend to exhibit an Anderson local¬ 
ization there. The enhanced flat feature of S{q,uj) ob¬ 
served in the random singlet state at A > Ac is cer¬ 
tainly consistent with such a picture. In this view, the 
randomness-irrelevant QSL state is an extended RVB- 
like state, while the randomness-relevant QSL state is an 
Anderson-localized state of singlets. Whether this naive 
picture captures correct physics behind the two types of 
QSL-like states needs further clarification. 

In contrast to the kagome case, the neutron-scattering 
data are not presently availabe for triangular organic 
salts due to the oganic nature of the material. An¬ 
other route to the random-singlet state of S' = 1/2 
triangular AFs might be an insulating mixed crystal 
Cs 2 Cu(Bri_ 2 ,Cla,) 4 . Experiments suggest that this com¬ 
pound might exhibit a non-magnetic ground state in the 
range x > 0.17—, at which the system possesses a consid¬ 
erable amount of randomness associated with the random 
arrangement of Cl“ and Br“. Interestingly, in this QSL- 
like regime, the gapless behavior including the T-linear 
low-temperature specific heat is observed^. It might be 
interesting to perform neutron-scattering measurements 
on this compound in its QSL regime to make a compari¬ 
son with our present data. 

The naive picture of the random singlet state, which 
is observed in common with the triangular and the 
kagome models, might be that tightly bound spin singlets 
are preferentially formed at stronger bonds, leaving 
loosely bound spin singlets or nearly free spins at weaker 
Jij bonds. Of course, such simple assignment of singlets 
to randomly distributed Jij bonds immediately meets 
contradiction or ‘frustration’, revealing that the singlet 
formation in the spatially random environment is not a 
trivial matter at all. A subtle balance between the poten¬ 
tial energy due to the local energy gain arising from the 
nonuniform Jy and the kinetic energy arising from the 
resonance of the local singlet states should determine the 
true ground state. In fact, in both cases of the triangular 
and the kagome models, we have observed that, although 
there generally exists a tendency that the strong singlet 
with larger negative • Sj-value tends to be formed at 
strong bonds with large -value, this tendency is quite 
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often violated in the sense that strong local singlet is 
sometimes formed at weaker bonds, or the singlet forma¬ 
tion remains weak even at stronger bonds. 

In summary, we investigated the nature of spin cor¬ 
relations in the randomness-induced QSL-like state, the 
random-singlet state, of the random-bond 5 = 1/2 AF 
Heisneberg model on the triangular and the kagome 
lattices by computing their static and dynamical spin 
structure factors by means of the exact diagonalization 
method. Gapless behaviors accompnied by the broad dis¬ 
tribution extending to higher energy, dependent on the 
wavevector only weakly, is observed in the dynamical spin 
structure factor 5(q, w) in the random-singlet states of 
the both models in common. Especially in the kagome 
case with strong randomness, 5(q, w) hardly dependeds 
on the wavevector q and exhibits an almost flat ditribu- 
tion in a wide range of uj, accompanied by the w = 0 peak. 
These features agree with the recent inelastic-neutron 
scattering data on a single-crystal hermertsmithite semi- 
quantitatively, suggesting that the QSL state observed in 
herbertsmithite is indeed the random-induced QSL state, 
i.e., the random-singlet state. 


ACKNOWLEDGEMENT 

The authors are thankful to Prof. K. Kanoda, Prof. 
H. Tanaka, Prof. I. Ono, Prof. C. Hotta, Prof. T. Sakai 
and Dr. K. Hiji for useful discussion. Prof. H. Tanaka 
pointed the possible importance of the Jahn-Teller distor¬ 
tion assiciated with the Cu^’*' in the triangular layer in as¬ 
signing the quenched randomness in the exchange inter¬ 
action of herbertsmithite. The ED calculation was per¬ 
formed by use of TITPACK Ver. 2. The QMC calcula¬ 
tion was carried out using the ALPS QMC application^^. 
They are thankful to ISSP, the University of Tokyo for 
providing us with CPU time. This study is supported 
by a Grant-in-Aid for Scientific Research No. 25247064. 
One of the authors (T.S.) acknowledges the financial sup¬ 
port from the Motizuki Eund of Yukawa Memorial Eoun- 
dation. 


Appendix A: 

In this appendix, as a typical model of unfrustrated 
5=1/2 random systems, we deal with the random-bond 
5=1/2 AE Heisenberg mdel on the square lattice. The 
Hamiltonian is given by eq. (1), the random interaction 
Jij obeying the uniform distribution characterized by the 
randomness parameter A, commonly with the triangu¬ 
lar and the kagome cases. The ED study on the bond¬ 
randomness effect in the unfrustrated square model have 
already been made in Refi^. It was shown there that, in 
the square-lattice case, the subalttice magnetization sur¬ 
vived against the bond randomness. In this appendix, we 
show not only the ground-state sublattice magnetization 
but also the temperature dependence of the sublattice 


magnetization, together with the static and the dynami¬ 
cal ground-state spin structure factors, in order to make 
full comparison with the corresponding data for the frus¬ 
trated systems. 

The lattice size is A = 8,10,16,18, 20, 24 and 32 
with periodic boundary conditions in all directions. The 
sample average is performed for 300 {N = 8,10), 200 
(A = 16,18), 100 (A = 20,24), and 5 (A = 32) inde¬ 
pendent bond realizations in calculating the sublattice 
magnetization, and for 100 (A < 32) and 20 (A = 32) 
independent bond realizations in calculating the static 
and the dynamical spin structure factors. 



1AfN 

Figure 10. (Color online) The rescaled T = 0 squared sublat¬ 
tice magnetization per spin, m^, of the square-lattice Heisen¬ 
berg antiferromagnet plotted versus 1 / y/N (A the lattice size) 
for several values of the randomness A. 



Figure 11. (Color online) The temperature dependence of the 
rescaled sublattice magnetization per spin, iris, of the square- 
lattice Heisenberg antiferromagnet for the size A = 36, and 
for the randomness A = 0 and 1.0. 

The size dependence of the T = 0 squared sublattice 
magnetization per spin, m^, is shown in Fig. 10 for vari¬ 
ous randomness A. The sublattice magnetization here is 
the one associated with the two-sublattice AF order. As 
can clearly be seen from the figure, rus remains nonzero 
even in the the rmodynamic limit for all values of A, 
indicating that the AF LRO persists up to the maxi¬ 
mal randomness of A = 1 as reported in ref. [82]. This 
is in sharp contrast to the cases of the triangular and 
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Figure 12. (Color online) Intensity plots of the T — 0 static 
spin structure factor S(q) of the square-lattice Heisenberg 
antiferromagnet in the wavevector {q^, qy)-plane for the ran¬ 
domness A = 0 (a), and 1 (b). The solid black line represents 
the zone boundary of 1st BZ of the square lattice. 



Figure 13. (Color online) The ai-dependence of the T = 0 
dynamical spin structure factor S(q,uj) of the square-lattice 
antiferromagnet taken at the wavevector q = (tTjTt) for the 
randomness A = 0 (a), and 1 (b). 


the kagome models where the AF LRO gives way to the 
random-singlet state for sufficiently strong randomness. 
Hence, not only the randomness but also the frustration 
plays a significant role in stabilizing the random-singlet 
state. In other words, all three elements, i.e., the strong 
quantum fluctuation, the frustration, and the quenched 
randomness, conspire to realize the present QSL state, 
the random-singlet state. 

In Fig.ll, we show the temperature dependence of the 
rescaled sublattice magnetization per spin, rhs, for the 
regular A = 0 and for the maximally random A = I 
cases. The lattice size is N = 36. The calculation here 
is made by use of the quantum Monte Carlo method^, 
which is possible due to the absence of frustration in the 
square-lattice model. In both the regular and the max¬ 
imally random cases, rhs increases monotonically with 
decreasing the temperature down to the lowest temper¬ 
ature studied, as can be seen from Fig.ll. In particu¬ 
lar, in contrast to the maximally random A = 1 case of 
the triangular model and the general case of the kagome 
model, the crossover behavior associated with a decrease 
of rhs in the lower temperature regime is not observed 
in the square mdoel even in the maximally random case. 
This suggests that the crossover behavior observed in the 
triangular and the kagome models might be a character¬ 
istic of the formation of the QSL-like state including the 
random-singlet state. 


The T = 0 static spin structure factor S (q) is shown in 
Fig. 12 as an intensity plot in the wavevector plane. The 
lattice size is iV = 24. The computed S'(q) exhibit rather 
sharp peaks at q = (±7r,±7r) corresponding to the AF 
LRO of the model both for A = 0 and A = 1. In fact, 
overall ^(q) looks quie similar between the regular case 
of A = 1 (a) and the maximal random case of A = 1 (b). 

The w-dependence of the dynamical spin structure fac¬ 
tor S{q,uj) at a wavevector q = (tt, tt) corresponding to 
the AF order is shown in Fig. 13 for the randomness A = 0 
(a), and for A = 1 (b). The dominant peak is a single 
magnon excitation, the peak location tending to zero to¬ 
ward the thermodynamic limit as shown in the insets. 
As in the static case, 5'(q, w) looks qualitatively similar 
between the regular and the maximally random cases, 
while the magnon peak becomes broader in the random 
case. 


Appendix B: 

In this appendix, we show the lattice shapes used in our 
ED calculation for various lattice sizes N. The cases of 
the triangular lattice, the kagome lattice and the square 
lattice are given in Figs.14, 15 and 16, respectively. In 
all cases, periodic boundary conditions are applied in all 
directions. 


A/= 18{ty pe2) A/=24{ty pe2) 



Figure 14. (Color online) The lattice shapes used in the exact 
diagonalization calculation of the triangular-lattice model for 
various N. Periodic boundary conditions are applied in all 
directions. 
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A/=12 A/=18(type2) 

-■o 




Figure 15. (Color online) The lattice shapes used in the exact 
diagonalization calculation of the kagome-lattice model for 
various N. Periodic boundary conditions are applied in all 
directions. 


N=24 N=8 



Figure 16. (Color online) The lattice shapes used in the ex¬ 
act diagonalization calculation of the sqaure-lattice model for 
various N. Periodic boundary conditions are applied in all 
directions. 
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